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TRANSPORT PROPERTIES AT HIGH TEMPERATURES OF 


C0 2 -N 2 -02-Ar GAS MIXTURES FOR PLANETARY 
ENTRY APPLICATIONS 

By Jerry S. Lee* and Percy J. Bobbitt 
Langley Research Center 

SUMMARY 

The equations and assumptions employed to compute the viscosity, thermal conduc- 
tivity, and Prandtl number for high-temperature gas mixtures which, at low tempera- 
tures, may be composed of arbitrary percentages of CO2, N2, O2, and Ar are presented. 
Ionization phenomena are not considered; hence, the results are limited to pressures and 
temperatures where dissociation reactions dominate. Numerical results are obtained for 
three different "atmospheric" gas mixtures: air, 100 percent C0 2 , and a mixture of 
43 percent CO2, 25 percent N2, and 32 percent Ar (by vol.). Comparisons of the latter 
two mixtures with air and, also, of air and CO2 with the results of other investigations 
are given. In addition, transport properties computed from the full equations of the first 
approximation of the Chapman-Enskog rigorous kinetic theory are compared to those 
computed from first and second approximations to these equations. 

INTRODUCTION 

The Martian atmosphere, as well as other planetary atmospheres, is believed to 
consist primarily of carbon dioxide, nitrogen, and argon. In order to perform meaningful 
boundary layer and wake analyses for entry spacecraft, reasonably accurate values of the 
high-temperature thermodynamic and transport properties of these atmospheres are 
required. Thermodynamic properties of several NASA Mars model atmospheres have 
already been calculated (refs. 1 to 4). In this report, attention is focused on the calcula- 
tion of viscosity, thermal conductivity, and Prandtl number of dilute, ideal gas mixtures 
in a state of local chemical equilibrium. The mixtures consist of a maximum of nine 
chemical species: CO2, O2, N2, NO, CO, O, N, C, and Ar. A specific mixture need not 
contain all nine of these species. 

There are several approaches to the computation of transport properties of multi- 
component mixtures. The first approach, such as presented by Kennard (ref. 5), makes 
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use of simple kinetic theory based on mean free path concepts. The second approach 
makes use of the more precise Chapman-Enskog theory (ref. 6) which was extended to 
multicomponent mixtures by Hirschfelder, Curtiss, and Bird (ref. 7) and which was modi- 
fied to account for multiple atomic interactions by Mason, Vanderslice, and Yos (ref. 8). 
The third approach makes use of the kinetic theory of polyatomic gas mixtures developed 
by Mon chick, Yun, and Mason (ref. 9). 

The first two approaches apply rigorously only to dilute, nonreactive, monatomic 
gas mixtures, and the third applies to dilute, nonreactive, polyatomic gas mixtures. 

Thus, in a strict sense, none of these theories can be used for the mixtures under con- 
sideration without modification. Of the three approaches, the third is the most applicable. 
However, utilization of this theory requires extensive knowledge of rotational and vibra- 
tional relaxation times which, unfortunately, are not currently available for all of the 
polyatomic species. In this work, a slight modification of the second monatomic gas 
theory, referred to as the rigorous kinetic theory, is used. 

The principal modification consists of the introduction of an expression (ref. 10) to 
account for the effects of internal degrees of freedom on thermal conductivity. For vis- 
cosity and the translational portion of the thermal conductivity, the equations resulting 
from the first approximation (i.e., the lowest Sonine approximation leading to nonzero 
results for the transport properties) of first-order Chapman-Enskog perturbation theory 
(ref. 7) are used directly. In other words, internal degrees of freedom and chemical 
reactions are assumed to have negligible perturbing effects on these equations. The 
results of Butler and Brokaw (ref. 11) are used to compute a reactive contribution to the 
thermal conductivity. It is also assumed that the collision integrals involved in all these 
equations, referred to herein as the full kinetic -theory equations, can be evaluated from 
orientation-averaged interaction potentials. 

Although the basic equations for calculating transport properties are well known 
(refs. 7 and 10 to 14), the collision integrals needed in these equations have not yet been 
determined for all pairs of species present in the atmospheric models. Fortunately, 
many of the interactions of importance also occur in air, and these have received con- 
siderable study (refs. 15 to 19). Herein, the tables of Yun and Mason (ref. 19), containing 
the most important air-collision integrals, are used. 

The procedure followed in obtaining collision integrals for the remaining interac- 
tions — namely those associated with CC>2, CO, C, and Ar — is a modification of the tech- 
nique employed by Amdur and Mason (ref. 20). For each interaction, appropriate low- 
temperature and high-temperature potentials are obtained from experimental data and/or 
approximate calculations. Then, the low-temperature and high-temperature collision 
integrals determined by these potentials are plotted as a function of temperature, and the 
two branches are joined by a smooth curve. Although this curve is not uniquely defined, 
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the technique should be markedly superior to using collision integrals generated from 
low -temperature potentials throughout the entire temperature range. 

Numerical results for viscosity, thermal conductivity, and Prandtl number are 
obtained for three different atmospheres: air, 100 percent CO 2 , and a mixture (referred 
to herein as the Mars model atmosphere) of 43 percent CO 2 , 25 percent N 2 , and 32 per- 
cent Ar (by vol.). Comparisons of the latter two atmospheres with air and of the present 
air results with those of Hansen (ref. 21) and Yos (ref. 22) are shown. The CO 2 results 
are compared with those of Thomas (ref. 23) and Van Tassell (ref. 24). Transport prop- 
erties obtained from the full kinetic -theory equations are compared with those obtained 
from first and second approximations to these equations. 

SYMBOLS 


A parameter in exponential repulsive potential (eq. (6)), electron volts 

Ajj collision integral ratio 

a parameter in exponential repulsive potential (eq. (6)), 1/angstrom 

collision integral ratio 
b impact parameter 

^l,ij’ < "2,ij ,< "3,ij’l constant defined by equation (BIO), (Bll), (B12), (B24), and (B25), 
c 4,ij» c 5,ij / respectively 


mixture specific heat at constant pressure (per mole of undissociated 
mixture) 

c n 

nondimensional mixture specific heat at constant pressure, — L 

R 




molar specific heat at constant pressure of species i 


'PA 


Cn i 

nondimensional specific heat at constant pressure of species i, 

R 


'ij 


binary diffusion coefficient of species i and j 


parameter in point center of repulsive force potential (eq. (5)) 
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d O’ d N 


E ik 



g 


AHj 


H ij 



hi 


k 





Mi 

M u 


n a 


distance from center of mass of molecule to O or N atom 

element in reactive conductivity determinants (eqs. (B30) and (B31)) 

nondimensional element in reactive conductivity determinants (eqs. (41) 
and (42)), ^ E ik 

initial relative speed of colliding particles 
enthalpy of reaction for reaction i 

AH t 

nondimensional enthalpy of reaction for reaction i, 

RT 

element in viscosity determinants (eqs. (B2) and (B3)) 

nondimensional element in viscosity determinants (eqs. (25) and (26)), r? 0 Hj,j 
molar enthalpy of species i 

h* 

nondimensional molar enthalpy of species i, — 

RT 

indices 

Boltzmann's constant 

elements in translational conductivity determinants (eqs. (B17) and (B18)) 

nondimensional elements in translational conductivity (eqs. (31) and (32)), 
^o E ij 

molecular weight of species i 
molecular weight of undissociated mixture 
mass per particle, grams 
Avogadro's number 
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n 


Np r 

N Pr,I 

N Pr,n 

n 

n ij 

P 

R 


T 

rpt 

t 

[X]i 

x i 

a 

a i 

“ij 

P 


Prandtl number 

first approximation to Prandtl number 
second approximation to Prandtl number 

parameter in point center of repulsive force potential (eq. (5)) 

stoichiometric coefficient of species j in reaction i 

pressure, atmospheres 

universal gas constant, ergs per °K-mole 

distance between colliding particles, angstroms 

value of r at potential energy minimum 

distance of closest approach of particles 

absolute temperature, °K 

T 

nondimensional absolute temperature, -~ 

e/k 

number of independent reactions 

symbolic representation of chemical symbol for species i 
mole fraction of species i 

parameter controlling steepness of repulsive part of exp-6 potential 

exp- 6 parameter for binary collision of particles of species i 

exp-6 parameter for binary collisions of particles of species i and j 

temperature parameter used in finding collision integrals for exponential 
repulsive potential (eq. (22)) 
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mill i 


ij 


A? 


ij 




V 


V 


''I 

* 

% 


r h 


% 

x 

X* 






nondimensional relative velocity, 



parameter in reactive conductivity equations 


nondimensional parameter in reactive conductivity equations, 



depth of potential well, electron volts 

depth of potential well for interaction of two particles of species i, 
electron volts 

depth of potential well for interaction of two particles of different species i 
and j , electron volts 

viscosity 

nondimensional viscosity, — 

Vo 

first approximation to viscosity 
second approximation to viscosity 
viscosity of pure species i 

reference viscosity, 3.20295 x 10“^ ff, grams per centimeter-second 

total thermal conductivity 

nondimensional total thermal conductivity, — 

^o 

first approximation to total thermal conductivity 
second approximation to total thermal conductivity 
frozen thermal conductivity, X^ + Xj n { 
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first approximation to frozen thermal conductivity 


second approximation to frozen thermal conductivity 


nondimensional frozen thermal conductivity, 




conductivity of pure species i 

internal contribution to total thermal conductivity 

X* rij" 

nondimensional internal conductivity, —— 

reference thermal conductivity, 4.16118 x 10 3 \Jf, 
ergs per centimeter-°K-second 

reactive contribution to total thermal conductivity 

A r 

nondimensional reactive conductivity, — 

^o 

first approximation to reactive conductivity 
translational contribution to total thermal conductivity 

nondimensional translational conductivity, — 

^o 

first approximation to translational conductivity 
second approximation to translational conductivity 
m-m- 

reduced mass, i — 1— 

mi + mj 

total number of species considered 


collision diameter, angstroms 



collision diameter for binary collision of species i 


Ojj collision diameter for binary collision of species i and j 

<fi interaction potential, electron volts 

X angle of deflection 

collision integral 


n(l,s) 

nondimensional collision integral, — — 

_(Z,s) 

“rigid sphere 

modified collision integral, o2f2^’ s ) , angstroms^ 


ANALYSIS 

Thermodynamic Properties 

The thermodynamic properties required in calculating mixture transport properties 
are the mole fractions, the mixture constant-pressure specific heat, the species constant- 
pressure specific heats, and the species enthalpies. For each mixture, these properties 
were calculated by using the computer program developed by Allison (ref. 3) and Newman 
and Allison (ref. 4). Figures 1,2, and 3 show the dependence of the mole fractions on 
temperature for pressures of 0.001, 0.01, 0.1, 1.0, and 10 atm for air, CO 2 , and the Mars 
model atmosphere (43%CC>2, 25% N2, 32% Ar (by vol.)), respectively. Figures 4, 5, and 6 
show the variation of the nondimensional specific heat of these mixtures with temperature 
for the same pressures. The nondimensional species specific heats and species enthal- 
pies are shown in figures 7 and 8, respectively, and apply for all three mixtures. 

Inspection of figures 1 to 3 reveals that electrons and ionized atoms become impor- 
tant at the higher temperatures, especially at low pressure. Since ionization is not 
included in subsequent transport-property calculations, it is necessary to specify regions 
of validity for these properties. More is said about this in the section "Results and 
Discussion." 


Potentials 

As noted in the Introduction, transport properties depend upon the collision integrals 
for each pair of chemical species in the mixture. These collision integrals, in turn, 
depend upon the interaction potentials for the possible binary encounters. It is therefore 
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appropriate at this point to enumerate all possible binary interactions for a mixture of 
CO 2 , O 2 , N 2 , NO, CO, O, N, C, and Ar. These interactions are given in table 1. 

With reference to table 1, parentheses around a pair designate that its interaction 
potential has been considered previously and that the corresponding collision integrals 
can be obtained from the tabular values of Yun and Mason (ref. 19). Brackets around a 
pair mean that an accurate interaction potential for the pair is not critical because one or 
both of the mole fractions associated with the pair are always small. Table 2 summa- 
rizes the potentials used in the computations. The atom-atom potentials listed in this 
table apply to the interaction of atoms in their ground electronic states. For the limited 
temperature range considered herein, such potentials are adequate. 

In the following paragraphs, some general comments are made about the interaction 
potentials. Then attention is focused on specific cases with special emphasis placed on 
pointing out the difficulties encountered and the assumptions invoked. 

The potentials used at low temperature are the Lennard- Jones 6-12 potential, 



and the modified Buckingham exp- 6 potential 



a 


( 1 ) 


( 2 ) 


Potential parameters (i.e., e, a, a, and rm) for many of the like-like interactions 
were obtained from references 7 and 25 for the 6-12 and exp-6 potentials, respectively. 

In the absence of experimental data, potential parameters for interactions of unlike mol- 
ecules were generated from those for like molecules by using empirical combining rules. 
For the 6-12 potential. 


°ij = 


°i + °j 


(3a) 



(3b) 


and for the exp- 6 potential, 


m ,ij 


+ r 


m,j 


(4a) 


2 


9 



1 




(4c) 


These rules hold well for the interaction of molecules i and j provided these mole- 
cules have essentially spherically symmetric force fields. The sources from which 
potential parameters were obtained are given in table 2. 

At high temperatures the repulsive portion of the potential primarily determines 
the collision integrals (ref. 7). Under these conditions, two potential forms were used: 
the point center of repulsive force potential 



and the exponential repulsive potential 


(5) 


0 = Ae" ar 


( 6 ) 


where d, n, A, and a are constants for a given interaction. Potentials derived from 
molecular beam scattering experiments are normally cast in the first form, although the 
second can be used. Results of semiempirical quantum mechanical calculations are 
usually presented in the second form which has some theoretical basis. According to 
Amdur and Mason (ref. 20), equation (6) usually remains valid over a wider range of r 
than equation (5). For both potentials, experiment seems to confirm a geometric mean 
rule for finding potentials for unlike atoms from those for like atoms (ref. 20). That is, 
if 0(D-D) denotes the interaction potential for two D atoms and if 0( E-E) denotes the 
interaction potential for two E atoms, then the interaction potential for atoms D and E is 
given by (refs. 20, 26, 27, and 28) 


0(D-E) = 1/0 (D-D) 0(E-E) 


(7) 


It should be emphasized that this rule does not apply to atom-molecule or molecule- 
molecule interactions. 

Many of the high-temperature atom -molecule and molecule-molecule interaction 
potentials listed in table 2 were obtained by averaging appropriate atom-atom potentials 
over all orientations as was done in references 19 and 28 for binary interactions of atoms 
with diatomic molecules and diatomic molecules with diatomic molecules. Appendix A 
extends this technique to include interactions with the triatomic CO 2 molecule. Specific 
cases are treated in the following sections. 
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Carbon monoxide interactions .- Figures 2 and 3 show that CO is a major constitu- 
ent in the CO 2 and C 02 -N 2 -Ar mixtures throughout most of the temperature range con- 
sidered. Consequently, reliable transport property computations hinge strongly on knowl- 
edge of potentials for the various CO-molecule and CO-atom interactions. Yet there 
exists only a meager amount of information on these potentials, especially at high tem- 
perature. Specifically, molecular beam scattering experiments on the Ar-CO system 
yield (ref. 29) 


4>(Ar-CO) = eV (2.09 A = r = 2.68 A) (8) 

r 6.99 

To the authors* knowledge, no high-temperature potentials are presently available for the 
remaining CO interactions. 

Because of this lack of data, it was necessary to resort to reasonable approxima- 
tions. The logic used proceeds as follows. Since CO and N 2 have very similar nuclear 
and electronic structures and since these factors play a dominant role in determining 
intermolecular potentials, it appears plausible to assume that CO and N 2 behave similarly 
during molecular encounters. An indication of the validity of this assumption is provided 
in figure 9 where the high-temperature scattering potential for Ar-N 2 * (ref. 29), 

0(Ar-N2) = eV (2.04 A = r = 2.53 A) (9) 

r 7.0o 

is compared to the Ar-CO potential given by equation (8). The percent difference in 
energy varies from about 3 percent at the smallest common value of r to about 4 per- 
cent at the highest common value of r. From the standpoint of transport-property com- 
putations, the comparison is extremely favorable because the potentials enter these com- 
putations solely through collision integrals, and these quantities are relatively insensitive 
to the magnitude of the potentials. Figures 10 and 11, which compare the two most 
important collision integrals for the potentials of equations (8) and (9), clearly demon- 
strate this point. The percent difference in corresponding collision integrals never 
exceeds 3.1 percent. Figures 12 and 13 compare low-temperature collision integrals for 
CO-CO and N 2 -N 2 interactions. The Lennard- Jones 6-12 potential with parameters 
derived from viscosity measurements (ref. 7) was used for each interaction. Again the 
agreement is remarkable. 

"''The high-temperature potential listed in table 2 for Ar-N2 differs from equation (9). 
Equation (9) is valid for a higher temperature range than considered in the transport- 
property computations. It is used here simply because its range of validity is very 
nearly the same as that of equation (8). Were this not true, the comparisons would 
be meaningless. 
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Based on these comparisons, it was assumed that the potential for each CO inter- 
action could be replaced by the potential for the analogous N 2 interaction; for example, 
CO-O 2 was treated as N 2 -O 2 . This assumption is incorporated in table 2. 

Argon interac tions.- With the exception of the Ar-C02, Ar-C, Ar-O, and Ar-NO 
interactions, table 2 provides adequate explanation of the methods employed to formulate 
interaction potentials for the various argon interactions. Potentials for Ar-C02 and 
Ar-C interactions are considered in later sections dealing with carbon dioxide and carbon 
interactions, respectively. This section deals with high-temperature potentials for the 
Ar-O and Ar-NO interactions. 

The high-temperature Ar-O potential was calculated from the geometric mean rule 
(eq. (7)) with 


0( Ar-Ar) = 3.23 x 10 4 e ' 4 - 464r eV 


( 10 ) 


and 


0(0-0) = 32.56 e -2 - 197r eV 


(ID 


where 0( Ar-Ar) was taken directly from reference 20 and 0(0-0) was obtained from 
a least-squares curve fit of Meador’s results (ref. 18). 

It should be noted that Vanderslice, Mason, and Maisch (ref. 17) also present a 
potential for the 0-0 interaction. When their result is used, 0(Ar-O) differs apprecia- 
bly from that given in table 2. However, the collision integrals corresponding to the two 
potentials differ very little. Meador’s result for 0(0-0) was chosen primarily because 
it leads to greater consistency with known low-temperature potentials when it is used in 
the construction of high-temperature potentials for more complex interactions, such as 
CO 2 -O 2 , by the orientation-averaging technique. The same situation exists for the N-N 
interaction, which arises in the computation of the Ar-N interaction. 


To calculate the high-temperature Ar-NO potential, the technique of averaging indi- 
vidual atom -atom interaction potentials over all molecular orientations was employed 
(ref. 19). This procedure yields 


a 0 e~ aAr ° r r \ 

)) = Y — (a Ar0 r + l) sinh (a Ar0 d 0 ) - 

4rO d O r 11 

A e -a ArN r 
+ ^ArJI |7 a ,_„ r + i\ 


a ArO^O cos h ^ a ArO^Cj) 

— 


( 12 ) 
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where A Ar0 , a Ar O’ A ArN’ and a ArN are ex P onential repulsive potential param- 
eters for the Ar-O and Ar-N potentials (eq. (6)); r is the distance from the center of 
mass of the nitric oxide molecule to the argon atom; and dQ and d^f are the distances 
from the center of mass of the nitric oxide molecule to its oxygen and nitrogen atoms, 
respectively. The Ar-O parameters can be obtained by inspection of <p( Ar-O) in 
table 2. The Ar-N parameters come from 

<M Ar-N) = 2025e -3 * 59r eV (13) 

Equation (13) was obtained from the geometric mean rule using the Ar-Ar potential pre- 
sented in equation (10) and Meador’s result for the N-N potential (ref. 18), 

<£(N-N) = 127e _2 - 716r eV (14) 


According to Herzberg (ref. 30), 


d Q = 0.53667 A 
d N = 0.61333 A 

In order to make use of available tables of collision integrals, a least-squares curve fit 
of equation (12) was performed to obtain the exponential repulsive potential presented in 
table 2. 

Carbon dioxide interactions.- All low-temperature potentials in table 2 for CO2 
interactions except those for CO2-CO2, CO2-NO, and CC>2-Ar are based on experimental 
diffusion coefficient determinations (ref. 31). Lack of such data for CO2-NO and CC>2-Ar 
interactions necessitated use of viscosity data and the combining rules given in equa- 
tions (3). 

The validity of using the low-temperature CO2-CO2 potential in table 2 in the calcu- 
lation of multicomponent mixture transport properties is open to question. Recent self- 
diffusion coefficient measurements for CO2 between 1103° K and 1944° K indicate that 
different sets of 6- 12 parameters are required to correlate viscosity and self-diffusion 
data (refs. 32 and 33). Use of the potential from viscosity measurements in the prediction 
of self -diffusion coefficients, and vice versa, can produce errors as high as 30 percent or 
40 percent, at temperatures around 2000° K. Quite likely the need for two distinct poten- 
tials to correlate experimental results arises from the fact that CO2 is a long linear mol- 
ecule with a nonspherical force field. Based on this evidence, one must conclude that the 
theory does not yield uniformly good results for CO2. 
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In an effort to decide whether the viscosity or self-diffusion potential was more 
appropriate for present purposes, the low-temperature viscosity and thermal conductivity 
of pure CC>2 were calculated by using both potentials. The viscosity potential led to good 
agreement with experimental viscosity and conductivity correlations (ref. 34) for CO2, 
and the diffusion potential led to poor agreement. Therefore, it was assumed that the 
viscosity potential should also be more appropriate for mixture computations. It should 
be noted, however, that the situation might be reversed in the computation of multicom- 
ponent diffusion-coefficient properties which are not considered herein. 

With the exception of the CO2-N and CO2-C interactions, the high-temperature 
C02-molecule and C02-atom potentials in table 2 were derived by averaging atomic 
interactions over all orientations. Appendix A contains all pertinent equations used in 
this procedure. Since these equations contain the C-O, C-C, C-N, and C-Ar interaction 
potentials and since these potentials are not currently available, it was necessary to 
introduce further approximations. The C-0 potential was approximated with the N-N 
potential; and in the last three interactions, the carbon atom was treated as a nitrogen 
atom. These approximations are not as critical as it might appear because the averaging 
process automatically emphasizes the interactions of the peripheral atoms. Finally, 
establishment of CO2-N2 and CO2-NO potentials requires knowledge of <f>( N-O). For 
both potentials, 


<p( N-O) = 69.17e" 2,556r eV (15) 

was used. This potential was found by curve-fitting Meador's results (ref. 18). After the 
introduction of these approximate potentials into the equations of appendix A, each 
orientation-averaged potential was replaced by a least-squares curve fit to obtain the 
results presented in table 2. 

Carbon interactions .- The assumption that carbon atoms can be replaced by nitrogen 
atoms was also utilized in obtaining potentials for the interaction of carbon atoms with the 
various atoms and molecules. This assumption is much more critical because carbon — 
other-atom potentials occur in every term of the orientation-averaged potentials. 


Collision Integrals 

As noted previously, interaction potentials enter transport-property computations 
through a set of collision integrals. For the interaction of particles i and j along a 
single potential (f > , the integrals are defined by the following formulas (ref. 7): 

dr 


X = n - 2b 


r\oo 


j. 


( 16 ) 




<P 


m 
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(17) 



cos^ 



In general, two ground state atoms can interact along any one of several potential 
energy curves. For this situation, equations (16) and (17) must be evaluated for each 
curve, and then all of the resulting collision integrals with the same superscripts are 
averaged with a priori assigned weighting factors to yield an effective set of collision 
integrals for the interaction (ref. 8).- These computations have been performed for the 
atom-atom interactions of importance in dissociating air (ref. 19), and fortunately, many 
of the same results also apply for the mixtures under consideration. 


The basic-mixture transport-property equations are usually expressed in terms of 
reduced collision integrals. In reference 7, which is probably the most widely used ref- 
erence for gaseous transport theory, the transport properties are formulated in terms of 
a set of quantities , defined as the ratio of equation (17) to the corresponding 

expression for the rigid sphere model; that is, 


qM* = 


sfo’ s ) 

rigid sphere 


(18) 


In the tabulations of the collision integrals of air, Yun and Mason (ref. 19) employed the 
slightly modified form, 


= a2s^’ s )* 


(19) 


where o denotes the collision diameter. Since these tabulations were used in the pres- 
ent work, this notation was adopted for all interactions. 

The following collision integrals and collision integral ratios are required to com- 
pute mixture viscosity and thermal conductivity: 


o}M> 

(20a) 

“ij 

(20b) 

n< 2 ’ 2) 

A* - lJ 
1] 

(20c) 

4 1 ’ 2 * . 4 1 ’ 3 ) 

ojW) 

“ij 

(20d) 
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In these equations, subscripts i and j designate quantities associated with the inter- 
action of species i and j. 

(1 (2 2) 

For a given pair of species, that is for fixed i and j, Opj ’ ’ and j ’ ' vary 
considerably with temperature. In addition, at a fixed temperature, both quantities differ 
significantly for different pairs of species. On the other hand, A*j and B*j vary little 
regardless of the temperature or of the pair of species considered. This fact is illus- 
trated in figures 14 and 15 where these quantities are plotted for the Lennard- Jones 6-12 
and modified Buckingham exp -6 potentials. Figure 16 shows the quantities for the expo- 
nential repulsive potential. The abscissa in figures 14 and 15 is defined as 

T* = — 7 - (21) 

e/k 

and that in figure 16 as 

^ 1 "( 4 ) (22) 

Generally, for the low-temperature potentials in table 2, the extreme right-hand portions 
of figures 14 and 15 apply. Furthermore, /3 normally ranges between 7 and 10. As can 
be seen from the graphs, little error should incur if A^ and B^j are arbitrarily 
assigned the constant values 


A*j = 1.18 

(23a) 

B ij = 1-15 

(23b) 


These values give slight emphasis to the exponential repulsive potential (fig. 16) because 
of its appropriateness at the higher temperatures of interest. At any rate, actual values 
should deviate from the assigned ones by less than 5 percent. These results were used 
in all transport-property calculations. 

h j\ (2 2) 

For Ojj ’ and ’ no simple approximations hold because of the significant 
variations mentioned; consequently, it was necessary to consider each individual inter- 
action. Collision integrals for interactions that occur in air were taken directly from 
reference 19. Collision integrals for the remaining interactions, given in table 3, were 
obtained by essentially the same method employed in references 19 and 20. The principal 
steps followed in this procedure were as follows. For each interaction, S2^’^ was 
computed for both the low- and high-temperature potentials and plotted as a function of 
temperature. This produced a curve consisting of two branches whose ranges of validity 
were estimated by using the technique of Hirschf elder and Eliason (ref. 35). The 
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were 


branches were then joined by a single smooth curve from which values of 

(2 2 ) 

read. This procedure was also used for ’ . Collision integrals for the Lennard- 

Jones 6-12 potential were obtained from reference 7; those for the modified Buckingham 
exp-6 potential from references 7 and 25; those for the point center of repulsive force 
potential from reference 36; and those for the exponential repulsive potential from 
reference 37. 

Although there is no way of checking the absolute accuracy of the results presented 
in table 3, it should be pointed out that the calculation procedure outlined is intrinsically 
more accurate for some interactions than others because of variation in the lengths of the 
"faired-in" curves. To illustrate this point, curves used for the Ar-C>2 and CO 2 -O 2 
interactions are shown in figures 17 to 19. Figure 20 gives for the CO 2 -O 2 

interaction computed from both the high- and low-temperature potentials. No faired 
curve is given, however, since the upper curve does not become valid until the tempera- 
ture exceeds the range considered. Clearly, the Ar-C >2 results should be more accurate 
than the CO 2 -O 2 results. Actually, the CO 2 -O 2 curves are typical of all curves encoun- 
tered for CO 2 interactions. This further bears out the point made earlier concerning the 
uncertainty associated with CO 2 interactions. 

Transport Properties 

In the preceding sections all requisite data for computing viscosity, thermal con- 
ductivity, and Prandtl number have been presented. This section summarizes convenient 
forms of the equations for these properties. In addition, possible causes of numerical 
inaccuracy are pointed out, and techniques for alleviating inaccuracy are given. 

Because of the complexity of the kinetic-theory equations for viscosity and thermal 
conductivity of multicomponent mixtures, investigators (refs. 38 and 39, for example) 
have developed a number of relatively simple approximations, which "normally" yield 
values within a few percent of those attained from the full equations. In this note, the full 
equations were employed for the following reasons: 

1. As pointed out by Brokaw (ref. 14), inherent deviations in individual properties 
introduced by the use of approximate formulas may be such that errors are magnified in 
forming such ratios as the Prandtl number. 

2. The systems under consideration are simple enough that the full equations can be 
solved very rapidly with the aid of a high-speed computer. 

For purposes of assessing the validity of using approximate equations and also of 
comparing the final results with those of other authors, first and' second approximations 
to the viscosity, translational thermal conductivity, and a first approximation to reactive 
conductivity are included in the following subsections. These equations as well as the 
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special forms of the full equations are derived from more familiar equations in 
appendix B. 

Viscosity.- It is shown in appendix B that the full kinetic -theory expression for the 
viscosity of a ^-component, monatomic gas mixture can be expressed in dimensionless 
forms as 


where 


V* = ' 


H*n 

TT* 

n 12 * * • 

TT^ 

U lv 

X 1 

H 12 

H* 

n 22 * * * 


x 2 

TT + 

n lv 

a 2v ’ * * 

Kv 

Xp 

X 1 

x 2 ... 

X V 

0 


* 1*20u ,2) 

h ** = t x: 


H ?j 


1 + 2 £ x i x j c l,ij n ij ’ 
j=l 


(24) 


(25) 




(i/j) (26) 


and Cj and C 2 are defined by equations (BIO) and (B 11). 

Also shown in appendix B are the first and second approximations to equation (24), 


v 2 


* YiL 

v i ~ A H?- 


and 


i=l H ii 


V 2 V V „ _ TT • 

u 


(27) 


i=l H ii i=l j=l 


(28) 
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Although equation (24) applies rigorously only to dilute monatomic gas mixtures, it 
is well known that the presence of internal degrees of freedom has very little influence 
on viscosity (ref. 7). Hence, this result was applied without modification to the mixtures 
of monatomic, diatomic, and triatomic species under consideration. 

In using equation (24), care must be exercised to obtain meaningful results. 
Inspections of equations (25) and (26) reveals that both H*j and H*j approach 0 as x^ 
approaches 0. Under these conditions, both the numerator and denominator in equa- 
tion (24) become small and the ratio approaches the indeterminant form 0/0. To apply 
the equation correctly, only those rows and columns associated with chemical species 
actually present in the mixture should be included. This means that the orders of the 
determinants may change with changes in temperature and pressure. In other words, 
when a given species disappears because of dissociation, its row and column must be 
omitted, and the remaining rows and columns must be shifted to fill the void. Likewise, 
when a new species appears, a new row and column must be included. However, the 
inclusion of new rows and columns can be avoided by initially formulating equation (24) 
so as to account for all possible species. Then one is concerned only with omitting 
appropriate rows and columns. This process can be actuated either by inspecting the 
magnitude of x^ or H- and was employed herein. 

In addition to providing a means of avoiding numerical difficulty, the deletion pro- 
cess inherently results in greater flexibility. Thus, a single general formulation 
including many species can easily be specialized to simpler mixtures with fewer species. 

Thermal conductivity . - As noted in the Introduction, the rigorous kinetic theory 
applies strictly to monatomic gases; thus, it does not account for the presence of internal 
degrees of freedom. Although these degrees of freedom do not significantly affect vis- 
cosity and diffusion, energy transition among the several degrees of freedom does mark- 
edly influence thermal conductivity. At present, however, available polyatomic mixture 
theories remain intractable so far as practical computations are concerned; therefore 
approximations must be invoked. Energy transport due to the combined effects of diffu- 
sion and chemical reactions can become extremely large under certain conditions; hence, 
it must be taken into account. 

All these effects can be included in approximate fashion by expressing the total con- 
ductivity as (ref. 10) 


X - Xf + Xj n t + X r (29) 

where X^ designates the translational contribution; Xi n t, the internal contribution; and 
Xr, the reactive contribution. A frozen conductivity is also used in the calculations and is 
defined as 
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Xf - Xt + A_int 


First and second approximations Af t and Xf^jj are determined by replacing Xt in the 
equation for Af with A t;I and A tj n, respectively. The following paragraphs set forth 
the equations used to compute each contribution. 

According to appendix B, the results of Muckenfuss and Curtiss (ref. 13) for the 
translational conductivity 2 can be expressed in nondimensional form as 
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L* = 

li 


'2xf^ii’ ) (Mi- ^ x i x j c 2,ij c 4,ij n ij 
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(31) 


L ij “ x i x j c 2,ij c 5,ij n ij ’ 


(i £ l) (32) 


The constants C 2 , C 4 , and C 5 are defined by equations (Bll), (B24), and (B25). 

The first- and second-order approximations to equation (30) are 

v x 2 

<1 - -4 l 4 - < 33 ) 

i=l L ii 


2 Two definitions are commonly used for translational conductivity. The present 
form applies when the energy flux is expressed in terms of the temperature gradient and 
the diffusion velocities. 


20 



and 



v 2 v v 


XiX .j L ij 


i=l L ii 


L — 1 L / T * T * 

i=l j=l MiLjj 

i£ 


(34) 


The discussion, given previously in connection with viscosity, concerning the 
appearance and disappearance of chemical species as temperature and/or pressure 
change also applies to the evaluation of equation (30). Equation (30) was initially formu- 
lated to include all nine species. For a mixture containing less than nine species, the 
correct form was obtained by first deleting appropriate rows and columns in the numer- 
ator and denominator determinants and then shifting the remaining rows and columns to 
fill the gaps. 

The internal contribution to the conductivity was computed from Hirschfelder's 
approximation (ref. 12) 


Mnt - 



together with his generalized Eucken expression (ref. 10), 


X i= X tl+ Mfei.5) 
i t,i t 1r 2/ 


(35) 


(36) 


For computational convenience, these equations were combined and the resulting expres- 
sion was written in nondimensional form as 



(37) 


The reactive contribution to the conductivity arises from the energy transport due 
to the combined effects of diffusion and chemical reactions. In regions of high tempera- 
ture, dissociation occurs and then the species diffuse down the resulting concentration 
gradient to cooler regions where recombination occurs. Simultaneously, there is a flux 
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of recombined species from lower to higher temperature. The net effect is a transport 
of chemical enthalpy from regions of high temperature to regions of lower temperature. 


The reactive contribution was computed from the original results of Butler and 
Brokaw (ref. 11).^ Their formulation requires that the chemical reactions needed to 
account for the species present in the mixture be expressed in a special way. From the 
total of v species present, t is chosen as independent. Then a system of t reac- 
tions is constructed such that each independent species occurs on the left-hand side of 
one, and only one, reaction and has -1 as its stoichiometric coefficient. Furthermore, 
the right-hand sides of the reactions contain dependent species only; that is, if njj 
denotes the stoichiometric coefficient of species j in the ith reaction and if (x]j 
denotes the chemical symbol of species j, the system of reactions takes the form 


1/ 

Mi= X uHj 

j=t+i 


(i = l,2,...,t) 


(38) 


Provided the reactions are expressed in this fashion, the nondimensional form of 
the reactive conductivity becomes (appendix B) 


4 = - 
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(40) 


^In reference 40, Brokaw gives a simpler formulation of this problem. However, 
the original formulation proved more convenient for present purposes. 
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The first approximation to x£ is given by 


. * Af 




i=l 
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(44) 
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This approximation corresponds to setting E^i ^ k) equal to 0 in equation (39). Thus, 
the approximation neglects coupling among the reactions in the sense that E*^ truly 
vanishes only when reactions i and k share no common dependent species. One could 
easily derive a second approximation to include coupling effects; however, since the 
approximations are used herein to aid in explaining deviations of present results from 
those of other investigators and since these investigators use equation (44), a second 
approximation to x£ was not considered. 

For the mixtures considered, accuracy problems were occasionally encountered in 
the evaluation of equation (39) when a single set of chemical reactions was used to 
describe the compositional changes throughout the entire ranges of temperature and pres- 
sure. This difficulty arises because some of the elements in both the numerator and 
denominator determinants become extremely large whenever one or more of the mole 
fractions approach 0. Under these circumstances, the compositional changes should be 
described by a different system of reactions. Thus, several different sets of reactions 
are needed to evaluate x£ over the entire ranges of temperature and pressure. The 
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technique used to automatically select appropriate reactions is discussed in the following 
paragraphs. 

When all nine species were present in significant quantity, the following set of 
reactions was used: 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

The number to the left of a reaction serves to identify the reaction and the independent 
component involved in the reaction; that is, CO 2 , O 2 , N 2 , NO, and C are designated as 
components 1, 2, 3, 4, and 5, respectively. The dependent components, CO, O, N, and Ar 
are labeled as 6, 7, 8, and 9, respectively. 

When one of the independent components vanishes while all the dependent compo- 
nents are present in significant amounts, it is obvious from equation (41) that the diagonal 
element corresponding to the reaction involving this component becomes unbounded. (For 
example, Eji — 00 as xj — 0^ Inspection of equations (41) and (42) reveals that all 
other elements assume correct values. Accordingly, the form of equation (39) appropri- 
ate to reactions (45), can be modified to give the correct form for this situation simply by 
deleting the row and column containing the unbounded element in each determinant and by 
shifting the remaining rows and columns. The process works just as well when several 
independent species vanish simultaneously. In this work, the deletion and shifting pro- 
cess was actuated by inspecting the magnitude of Ej^. 

When one or more of the species chosen as dependent vanish, the situation cannot be 
handled as easily. It is not sufficient to omit reactions in equations (45); rather, com- 
pletely new sets of reactions must be formulated. This is true because some species 
chosen as independent in equations (45) become dependent under these circumstances. 

The following cases proved sufficient for the three mixtures considered: 

Case I Equations (45) 

Case II C0 2 Z CO + O 

o 2 z 20 
no z 4 n 2 + o 

a 


co 2 ~ co + o' 

0 2 Z 20 
N 2 Z 2N ) 

no;n + o 
c:co-o 

-J 


(45) 
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Case HI 


co 2 ~co+| o 2 

NO ~ | N 2 +|o 2 

Case IV No reaction 


For each case, vanishing of an independent species was handled by the technique 
previously explained for case I. Selection of the appropriate case was based on inspec- 
tion of the mole fractions. 

Prandtl number.- The Prandtl number is defined by 


Np r = 


*7C p 

M U X 


(46) 


In terms of the variables used herein, this equation becomes 


N Pr = 


16 ^ C P 
25 M U A* 


(47) 


where Mu denotes the molar mass of the undissociated mixture in g/mole. The Prandtl 
number defined by equation (47) is frequently termed the ’’effective Prandtl number” to 
connote that both c* and X* include reactive contributions. Two approximations to 
the Prandtl number Np r j and Np r were computed for air and C0 2 • These approx- 
imations are given by 


N Pr,I = 


16 VP 

25 M u Xj 


(48) 


where 
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X I = X t,I + X int + X r,I 
X H ~ X t,H + X int + X r,I 
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RESULTS AND DISCUSSION 


In this section numerical results are given, in turn, for air, carbon dioxide, and the 
Mars model atmosphere (43% CO 2 , 25% N 2 , and 32% Ar (by vol.)). For each mixture, 
calculations based on the full equations are presented first. Then, for air and carbon 
dioxide, calculations based on the full equations are compared both with those based on 
approximate equations and with those of other investigators. Finally, the viscosity, total 
conductivity, and Prandtl number of air are compared to the respective properties of 
carbon dioxide and the Mars model atmosphere. 

Since ionization was not taken into account in the computations, the results at a 
given pressure are terminated at a cut-off temperature which was estimated by com- 
paring thermodynamic properties, Cp in particular, computed with and without ioniza- 
tion taken into account. Normally, the temperature cut-off determined in this manner 
corresponds to an electronic mole fraction between 5 x and 10“^. 

Transport Properties of Air 

The full-equation computations for viscosity, frozen conductivity, reactive conduc- 
tivity, total conductivity, and Prandtl number are shown in figures 21, 22, 23, 24, and 25, 
respectively. For a given pressure, the temperature dependence of each property can be 
divided roughly into three portions, corresponding to the reactions taking place. From 
1000° K to the region at which the curves start to diverge, the mixture remains essen- 
tially chemically inert. Divergence can be attributed largely to the dissociation of oxy- 
gen and to the formation and subsequent dissociation of nitric oxide. The large change in 
properties at still higher temperature arises primarily because of nitrogen dissociation. 

Figure 26 compares the viscosity from full-equation calculations at p = 1.0 atm 
with the first and second approximations and also with the calculations of Hansen (ref. 21) 
and Yos (ref. 22). The curve for the first approximation falls consistently below that 
for the full-equation calculation, reaching a maximum deviation of approximately 
11.5 percent at 6500° K. The second-approximation curve lies below but much closer to 
the full-equation curve; in fact, the deviations never exceed 1.5 percent. Yos' results lie 
slightly above the full-equation values at low temperatures and slightly below at high 
temperatures. Throughout most of the temperature range, Yos' values agree well with 
the second-approximation values. This agreement was expected because he used the 
same collision integrals for the major interactions and Brokaw's approximate equation 
(ref. 38) derived from a third-order approximation to the full viscosity equation. In con- 
trast, the agreement with Hansen's calculations is much worse, especially above 1500° K. 
At the higher temperatures, his results consistently deviate from the full-equation values 
by about 23 percent to 24 percent. The reason for this marked difference can be 



attributed to the fact that he uses simple kinetic theory and rough approximations for 
collision cross sections. 

As shown in figure 27, the same type situation exists for the frozen conductivity at 
p = 1.0 atm as for viscosity. The curve for the first approximation always lies below 
that for the full-equation calculation and has a maximum deviation of about 19 percent 
between 6000° K and 7000° K. The second -approximation curve is also low, but deviates 
by less than 6 percent. Yos’ results, which were calculated from Brokaw's formula 
derived from a third-order approximation, are in excellent agreement with the full- 
equation values. On the other hand, Hansen's results are again low with maximum devia- 
tion of about 29 percent at 6500° K. 

Figure 28 illustrates the danger inherent in neglecting coupling among the reac- 
tions, or equivalently, in using the first approximation, in the computation of the reactive 
contribution to the conductivity. Although the figure does not show it clearly, there is 
substantial disagreement in the full-equation and approximate values near 2000° K. This 
disagreement is relatively unimportant, however, because the reactive conductivity is a 
small portion of the total conductivity in this temperature range. On the other hand, the 
differences which arise from neglecting coupling due to NO dissociation between 4000° K 
and 6000° K are extremely significant because Ay is the principal part of the total con- 
ductivity in this temperature range. The differences are even larger than they appear at 
first glance. For instance, the percentage deviation of the first-approximation curve 
reaches a maximum of 61.6 percent at 4700° K. 

Also shown in the same figure are the calculations of Yos (ref. 22) and Hansen 
(ref. 21). Since Yos used the first approximation to the reactive conductivity and Yun and 
Mason's collision integrals, his result should, in principle, lie precisely on the dashed 
curve. To discern why the point at 5000° K does not coincide with the dashed curve, Yos' 
sources of thermodynamic information were investigated. Yos used mole fractions from 
his reference 43 for the temperature range from 1000° K to 5000° K and mole fractions 
from his reference 44 for temperatures above 5000° K. He also employed unpublished 
calculations made at AVCO Corporation for his Cp ^R and h^RT. It was foupd that 
his reference 43 values of the mole fraction for NO differed greatly from the values used 
herein, the largest deviation being 80 percent. As a check, the first approximation was 
recalculated using Yos* sources of mole fractions and the present Cp ^R and Iq/RT. 
Excellent agreement with Yos' values of were then obtained throughout the entire 
temperature range. Hansen's results agree well with the first-approximation values up 
to about 3000° K. His exceptionally low values between the two peaks are due primarily 
to his complete neglect of the reaction NO Z N + O. The lowness of Hansen's values in 
the vicinity of the second peak is thought to be a direct consequence of his method for 
obtaining collision cross sections because it has been shown that uncertainties associated 
with this method increase rapidly at high temperature (ref. 23). 
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The same general comments apply equally well to the total conductivity shown in 
figure 29. Yos' results should agree closely, but not exactly, with the second- 
approximation results, whereas Hansen’s results should be in better agreement with the 
first-approximation results. It is especially noteworthy that the second-approximation 
value deviates from the full-equation value by 29.5 percent at 4700° K. Thus, the pre- 
viously noted substantial divergence of the approximate reactive conductivity perpetrates 
into the approximate total conductivity in a somewhat weakened but still significant 
manner. 

Since viscosity and conductivity enter directly into the Prandtl number, it is obvious 
that the Prandtl number can be extremely sensitive to approximate calculations. This is 
borne out in figure 30, which shows the full-equation Prandtl number, the first- and 
second-approximation Prandtl numbers, and the results of Yos (ref. 22) and Hansen 
(ref. 21). The striking differences in the approximate and full-equation calculations are 
due primarily to the reactive -conductivity approximation. The largest deviations, occur- 
ring at 4700° K, are 25.4 percent and 23 percent for the first- and second-approximation 
values, respectively. Comparisons with the results of Yos and Hansen are complicated 
by the fact that their values reflect differences in thermodynamic data. Hansen’s values 
of Cp do not differ markedly from those used herein, and accordingly, most of the dif- 
ference shown in the figure can be attributed to the use of different theories and to his 
neglect of the NO reaction. In view of the rather large disagreement noted previously 
for ri* and A*, it is indeed surprising that Hansen’s Prandtl numbers agree with the 
full-equation values as well as they do. As noted previously, Yos employed two different 
sources for mixture composition and one of these sources gives NO mole fractions which 
differ markedly from the present ones. His calculation of the Prandtl number is incon- 
sistent with these sources in the following sense. Yos used the composition from his 
reference 43 for temperatures between 1000° K and 5000° K and that from his refer- 
ence 44 at higher temperatures, but he employed the c p ^R values from his reference 43 
throughout the whole temperature range. In order to obtain a useful comparison with his 
results, the second approximation to the Prandtl number was recomputed using Yos’ 
sources of mole fractions and c p ^R and the present c p ,i/ R and hi/RT. Figure 31 
shows that these new values are in good agreement with Yos’ calculations. It is evident 
that the differences between Yos’ Prandtl number and the second-approximation values 
are primarily due to the differences in the thermodynamic data. 

Transport Properties of Carbon Dioxide 

The results of the full-equation calculations of viscosity, frozen conductivity, reac- 
tive conductivity, total conductivity, and Prandtl number of carbon dioxide are presented 
in figures 32, 33, 34, 35, and 36, respectively. Just as with air, for a fixed pressure, the 
temperature variation of each of these properties can be split roughly into three regions 
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depending upon the reactions taking place. Up to the point at which the curves initially 
diverge, the gas consists of practically pure CO 2 . As the temperature increases beyond 
this point, CO 2 starts to dissociate into CO and O 2 , and with a small additional increase 
in temperature, O 2 also starts to dissociate. Thus, the initial divergence, and for some 
curves, the first peak, arises from the almost simultaneous dissociation of CO 2 and O 2 . 
After these reactions essentially go to completion, CO dissociation begins. This reaction 
is responsible for the rapid rise in viscosity and frozen conductivity and for the second 
peaks in reactive conductivity, total conductivity, and Prandtl number. 

For each of these properties, results of the full-equation calculations are compared 
with approximate computations and with the results of Thomas (ref. 23) and Van Tassell 
(ref. 24) in figures 37 to 41 for a pressure of 0.1 atm. Both Thomas and Van Tassell 
attained widely varying Prandtl numbers at this particular pressure; therefore, this pres- 
sure was chosen to ascertain why this variation was true. (Thomas also attained wide 
variation at lower pressures but Van Tassell did not present data for p < 0. 1 atm.) 

In figures 37 and 38, the approximations for viscosity and frozen conductivity com- 
pare with the full- equation calculations in much the same fashion as was found in the 
comparison for air. Van Tassell's viscosity calculations were expected to be in good 
agreement with either the exact values or the second approximation because he utilized 
Yos T method (ref. 22) which was previously found to be quite good for air. No data from 
Van Tassell is included in figure 38 because he does not present results for frozen con- 
ductivity. For both viscosity and frozen conductivity, Thomas' results lie considerably 
below even the first approximation; this too was expected since he utilized Hansen's 
approach (ref. 21). 

Figure 39, showing the reactive-conductivity comparisons, is particularly illumi- 
nating. In the region between 1500° K and 2700° K, the approximation values are much 
too low; in fact, the deviation from the solid curve approaches 100 percent near 1500° K. 
Between 2700° K and 3500° K, the approximation yields exceptionally high predictions, 
with a maximum deviation of 43 percent near 3000° K. No data are shown from Van 
Tassell because he does not provide sufficient data to allow determination of this con- 
tribution. With the exception of the point at 3000° K, Thomas' results agree well with 
the exact calculation up to 5000° K. The lowness of his values in the vicinity of the 
peaks is thought to be due to two possible causes: 

(1) As noted earlier, uncertainties involved in his method of finding collision cross 
sections become larger with increasing temperature. 

(2) Although he used the full equation for X r , his equation for does not agree 

with the original source (ref. 11). It is not clear whether this inaccurate equation was 
actually used in his computations. 
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Figure 40 shows clearly that the behavior of the reactive conductivity carries over 
directly to the total conductivity. Only the first-approximation curve is shown because 
the first- and second-approximation curves are practically indistinguishable. The data 
point of Thomas at 4000° K seems to be in error because his tables for A, Af, and A r 
are inconsistent at this temperature. 

Figure 41 clearly shows that Van Tassell’s high values for Prandtl number around 
2000° K are directly attributable to the use of approximate equations. In particular, the 
exhibited behavior is directly traceable to the reactive-conductivity approximation. The 
large disagreement of Thomas' values at 4000° K and 5000° K with the exact values is 
attributable, respectively, to a tabular inconsistency and to an extremely low specific 
heat. 


Transport Properties of the Mars Model Atmosphere 

The viscosity, frozen conductivity, reactive conductivity, total conductivity, and 
Prandtl number, as computed with the full equations, are presented in figures 42, 43, 44, 
45, and 46, respectively. Qualitatively, these properties exhibit the same type tempera- 
ture variation as previously observed for air and carbon dioxide. There are, of course, 
distinct quantitative differences and these are treated in detail in the next section. 

The initial divergence of the constant pressure curves given for these properties 
stems from the dissociation of CO 2 . With a slight increase in temperature, O 2 dissocia- 
tion also enters the picture. Thus, the lower temperature portions of the curves are 
influenced by the coupled dissociation of CO 2 and O 2 . As these reactions tend toward 
completion, CO and N2 dissociations become predominant’. It is these two reactions 
which are largely responsible for the steep increase in viscosity and frozen conductivity 
and for the second set of peaks in reactive conductivity, total conductivity and Prandtl 
number. Although argon does not react to any appreciable extent, it does affect the quan- 
titative character of the transport properties both through molecular interaction and 
through its influence on the equilibrium composition. 

Although no approximate calculations are given for this mixture, it is evident from 
the two preceding examples that use of approximate equations would lead to serious error 
in the total conductivity, and hence, the Prandtl number. This observation would be espe- 
cially true in the low temperature range where CO 2 and O 2 dissociate simultaneously. 

To the authors’ knowledge, transport properties for this mixture have not been 
computed previously; hence, no comparisons are possible. 
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Comparisons of the Transport Properties of Air, Carbon Dioxide, 
and the Mars Model Atmosphere 

In this section, the three mixtures are compared on the basis of the full-equation 
calculations of viscosity, total conductivity, and Prandtl number. For convenience, air 
is adopted as a standard and the properties of carbon dioxide and the Mars model atmo- 
sphere are discussed in terms of their deviation from the corresponding air properties. 
Attention is restricted to a pressure of 1.0 atm because similar results hold for the other 
pressures. 

Figure 47 compares the viscosities of the three mixtures. The viscosities of air 
and carbon dioxide are quite similar in magnitude, the viscosity of the latter never devi- 
ating from that of the former more than 6 percent. The viscosity of the Mars model 
atmosphere always lies above that for both air and carbon dioxide, but its deviation from 
that of the former is always less than 12 percent. There are, of course, distinctive fea- 
tures in the fine details of the temperature variations, such as the characteristic rises 
associated with the occurrence of chemical reactions. 

The curves for air and carbon dioxide are expected to be similar after appreciable 
dissociation of carbon dioxide, because the collision integrals and thermodynamic prop- 
erties of the remaining atoms and diatomic species are similar. The higher viscosity of 
the Mars model atmosphere is due mainly to the large amount of argon present in this 
mixture. 

Figure 48 illustrates the tremendous differences in the total conductivity of the 
mixtures. In the temperature region between 2000° K and 4000° K, the conductivity of 
CO 2 deviates from that of air by as much as 180 percent and between 5000° K and 6000° K 
by as much as 70 percent. For the Mars model atmosphere, deviations range from about 
95 percent near 2500° K to about 34 percent near 5500° K. 

As shown in figure 49, these large differences in total conductivity do not influence 
the Prandtl numbers as much as one might expect. The reason for this, of course, is the 
analogous behavior of the specific heat and total conductivity of the individual mixtures. 
The deviation of the Prandtl number of CO 2 from that for air ranges from about 12 per- 
cent to 20 percent. The maximum deviation of the Prandtl number of the Mars model 
atmosphere from that of air is about 10 percent. According to Horton and Zeh (ref. 41), 
such difference can be significant in problems such as blunt-body stagnation-point heat 
transfer. 

It is noteworthy that Prandtl number comparisons based on the second approxima- 
tion would differ greatly from those shown in figure 49. In particular, deviations of sev- 
eral hundred percent would be observed. 
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CONCLUSIONS 


On the basis of the foregoing analysis, the following conclusions are justified: 

1. Neglect of coupling terms in reactive-conductivity calculations for mixtures 
containing CO 2 can induce order of magnitude errors in total conductivity, and hence, 
large errors in Prandtl number. 

2. For the mixtures considered, second approximations give reasonable accuracy 
(§6 percent deviation from full-equation calculations) for viscosity and translational 
thermal conductivity. 

3. There exists a definite need for accurate interaction potentials for carbon-atom 
and carbon-molecule interactions. The accuracy of high-temperature transport proper- 
ties of mixtures containing carbon elements and/or compounds must be regarded as 
uncertain until such potentials become available. However, the present calculations are 
regarded as the best presently available. 

4. Accurate mixture specific heat is essential to accurate Prandtl number 
calculations. 

It should be emphasized that conclusions 1 and 2 are independent of conclusion 3 
because carbon — other -particle interaction potentials were used consistently in all 
calculations. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 1, 1969. 
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APPENDIX A 

ORIENTATION-AVERAGED POTENTIALS FOR C0 2 INTERACTIONS 

The purpose of this appendix is to develop high -temperature potentials for binary 
interactions of CO2 with atoms, diatomic molecules, and other C0 2 molecules. The 
basic idea is to construct the potential for two interacting particles from known potentials 
for interactions of atoms comprising the particles. This leads directly to an orientation- 
dependent particle potential which is then averaged over all possible orientations to yield 
an effective potential dependent only on the center-of-mass separation of the particles. 
Appropriate equations for atom — diatomic -molecule and diatomic -molecule — diatomic- 
molecule potentials are presented in reference 19. 

The procedure incorporates the following simplifying assumptions (ref. 19): 

(1) Each molecule can be treated as a rigid structure with ground-state dimensions. 

(2) Atomic interactions are independent and each atomic interaction potential has 
the form 

<t> = Ae-* 1 *- (Al) 


where R is the interatomic separation distance. 

(3) All molecular orientations are equally probable. 

Sketch 1 shows appropriate coordinate systems and nomenclature for the establish- 
ment of the interaction potential for the encounter of molecule AB with a symmetric, tri- 
atomic linear molecule DCD. For a particular configuration, assumption (2) implies that 

0(r,0i,a>i,0 2 ,w 2 ) = <£ ac + <t> BC + 0AD* + ^AD 2 + ^BDj + ^BD 2 ( A2 ) 


Xl X 2 
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where each term has the form given in equation (Al). Subscripts 1 and 2 are affixed 
to the D atoms because the distances from these atoms to atoms A and B generally differ, 
and hence, the potentials also differ; that is, in general, / ^AD2 ancI 

^BDj ^ ^BD 2 ‘ 

The orientation-averaged potential, <p(r) is related to $(r,0i,ct>i,0 2 > w 2) by 


1 p2ir /» 7T • p7T 

0(r) ^ \ \ \ \ 0 (r, 0 i,w 1 , 02 J w 2 ) sin sin 6» 2 d(9i d0 2 dw x du> 2 

(4^)^ J co 2 =0 J0 1 =O X 


Note that no Boltzmann weighting factor occurs in this equation; this is consistent with 
assumption (3). 

When equation (A2) is substituted into equation (A3), one finds that the six resulting 
integrals are quite complex and difficult to integrate. However, one can easily demon- 
strate that the integral of $ A( -. ^ as same f° rm as that for 4>j$c and also that the 
integrals of ^ , 4> ^ <P BDl ’ and ^BD 2 ^ ave the same form. Hence (p( r) can 

easily be found once the averaged potential for the A-C and A-Di interactions are known. 


. The averages of <p^ c and (p^ 
must obviously be independent of the par 
ticular choice of coordinate systems. 
Sketch 2 shows a new coordinate system 
which leads to marked simplifications in 


the integrals. Denoting the averaged val- 
ues of 0 AC and (p^ by ^> AC ^ and 



, respectively, one obtains 



<^Ac) = "f^I 0 e aAcRAC sin &i d6 1 

^>ad^ = J 0 e aAI)RAI) sin e l sin 0 2 d©i de 2 

where 

/ _ \l/2 

R ac = (r2 - 2rd A cos $i + d%) 


(A4) 

(A5) 

(A6) 
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( o 2 ^ ^/ 2 

R AD = ( R AC “ 2R AC d D cos 0 2 + d D/ 


(A7) 


Carrying out the indicated operations yields 


jAl 0 i 

($AQ) ~ ( a AC r + *) sinh ( a AC d A) “ a AC d A cosh ( a AC d A) 

' 7 at^d.r LI J 


-a AC r 


a AC d A r 


{^ad} “ 


A AD e 


- a AD r 


a AD d A d D r 


( a AD r + 2 ) Sinh ( a AD d A) sinh ( a AD d D) 


- a AD d A cosh ( a AD d A) sinh ( a AD d D) 

- a AD d D cosh ( a AD d D) sinh ( a AD d A)J 


(A8) 


(A9) 


Averages for the remaining atomic potentials in equation (A2) can now be obtained by 
changing the subscripts appropriately in equations (A8) and (A9). 

Use of these results leads directly to the following potentials for the CO2-NO, 
CO2-O2, and CO2-N2 interactions: 


A„„e' a NC 


<p(c0 2 -N0j " a 2 d r j( a NC r + *) sinh ( a NC d N) " a NC d N cosh ( a NC d N)j + ^ ~ j( a OC r + *) sinh ( a OC d o) ' A OC d O cosh ( a OC d C>)] 


a NO d N d D 


— [( a NO r + 2 ) sinh ( a NO d N) Sinh ( a NO d D) - a NO d N cosh (Wn) slnh (Wd) * a NO d D cosh (Wd) slnh 




OOOD 


r |( a oo r + 2 ) sinh ( a oo d o) sinh ( a oo d r) " a oo d o cosh ( a oo d o) sinh ( a oo d r) ~ a oo d D cosh ( a oo d p) sinh ( a oo d o)J (A 10) 


<p(cOz-o£j - — — [( a oc r + - 1 ) sinh ( a oc d cJ) ~ a oc d o cosh ( a oc d o)j + — “l( a OO r + 2 ) sinh ( a oo d o) sinh ( a OO d D) 

a nr d n r L L 


2 

a OC u O J 


D d O cosh ( a OO d c) slnh ( a OO d lJ " a CX) d D Cosh ( a OO d D) sinh ( a OO d C^| 


(All) 
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0(co 2 -n 2 ) 


2A NC e 


" a NC r 


a NC d N r 



l ) 


sinh 



a NC d N cosh 



4A e " aN ° r i 

i 4A NO e 

O 

( a NO r + 2 ) Sinh ( a NO d N) Sinh ( 

a NO d D) 


a NO d N d D r 1 

LI 

\ * 



" a NO d N cosh ( 

a NO d n) sinh 1 

a NO d o) 



" a NO d D cosh ( 

»NO d c) 8lnh | 

a NO d N)] 


(A 12) 


In each of these equations d D denotes the length of the C-O bond in the C0 2 molecule. 
The meaning of the other d symbols depends upon the particular interaction. For 
example, in equation (A10) d N denotes the distance of the N atom from the center of 
mass of the NO molecule whereas in equation (A12) d-^ denotes the distance of the 
N atom from the center of mass of the N 2 molecule. 

By similar analyses one can also establish the following potentials: 


0^CO 2 -CO 2 J — A£ C e 


- a CC r 


4 A oc e 0 0 
a OcV 


4A e " a OO r 

* A oo 


a OO d D r 


j( a OC r + *) slnh ( a OC d D) • a OC d D cosh ( a OC d D)] 
( a OO r + 2 ) Sinh ( a OO d D) Sinh ( a OO d D) 


- 2a OO d D C0Sh ( a OO d D) Sinh ( a OO d D) 


(A13) 


< p(c o 2 -n) = 


A " a NC r 2A NO e 

a nc + 


a NO r 


a NO d D r 


j( a NO r + *) Sinh ( a NO d D) " a NO d D C0Sh ( a NO d I^j 


(A14) 
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- A e " a °C r + 2A OO e 

- A OC e + Y 


a OCrD 


~ a OO T f- 

sHl' 


a oo r + 



a OO d D cosh 



(A 15) 


C °2-Ar) = A ArC e" aArCr 


2A ArO e 


' a ArO r 


a ArO d D r 


|( a ArO r + X ) sinh ( a ArO d D) “ a ArO d D cosh 



(A16) 
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ESTABLISHMENT OF EQUATIONS FOR VISCOSITY, TRANSLATIONAL 
CONDUCTIVITY, AND REACTIVE CONDUCTIVITY 

In this appendix, the expressions for viscosity and translational thermal conductiv- 
ity resulting from the first approximation of first-order Chapman-Enskog kinetic theory 
are cast into convenient forms for numerical computation. A convenient form for the 
reactive contribution to the thermal conductivity is also derived, and a first and second 
approximation for viscosity and translational conductivity and a first approximation for 
the reactive conductivity are developed. 

The viscosity of a I'-component gas mixture is given by (ref. 7) 


H 11 H 12 ' * * H ly X 1 

H 12 H 22 ‘ * * H 2y x 2 


®1j/ 

Xj X 2 ... x v 0 


where M denotes the determinant of the Hy elements defined by 


H _*i , V 2x i x .i RT L 3 M j * 
Hu " \ + 4 Mj + Mj ^l 1 + 5 Mi Al J 




(Bl) 


(B2) 


H 


i] 


2x i X j RT L 3 p*\ 
Mj + Mj pDjj ^ 5 7 


(i t j) (B3) 
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Here 77^, the viscosity of species i, is given by 

a/2 




Vi = 


5 / R 


16 K, 


cm. 


j/MjT 

“ii 


(B4) 




with expressed in units of A^. The binary diffusion coefficient of species i 

and j, Djj, is given by 


D ij= < 



\2 

10l6(^_ 

cm/ 


fgq + Mj)t 3 


vd}' 1 '! 2MiM j 


(B5) 


“i] 




where is also expressed in A^. 

It is convenient to arbitrarily define a reference viscosity tj 0 by 

r , r- .-n 

'3/R \* 2 X 


Vo = 


|8 K, 


loie 


g 


\cm -molef 


\ 1/2 

' iff 


(B6) 


This particular choice was made solely because it leads to convenient numerical values 
for the dimensionless viscosity rj*. Dividing equation (Bl) by equation (B6) and substi- 
tuting equations (B4) and (B5) into equations (B2) and (B3) give 


* = _7_ = 
Vo 


H *ll 

* 

h 12 

• • 

* 

x i 

h ! 2 

H 22 

* * 

• h 2- 

x 2 

• 

* 

• 

* 

• 

• 

* 

• 

* 

• 

• 

* 

• 

* 

• 

K lv 

H 2, 

• • 

• Ku 

x„ 

X 1 

x 2 

• * 


0 



H *i 




(B7) 
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where 


1.2J2< 2 ’\ 2 f -(1,1) 

" + 2 jli X,X,Cl ’ I,ail 

j£ 

(B8) 

H Ij'-2xiX j C 2iij 4 1 ’ 1) (l-f^) 

(i / j) (B9) 

Ci ij = (l +1— 1 A^C 2 ji 
1)13 l 5 VJ ^> 13 

(BIO) 

c Cs,i3 

2 ’ 13 Mi + Mj 

(Bll) 

2M i M j 

° 3 ’ 13 " Mi + Mj 

(B12) 


The reasons for defining C in this fashion are twofold. First, they can be computed 
once at the beginning of the calculations and used throughout the whole range of tempera- 
tures and pressures. ^Recall that A*j is taken as being constant.) Second, some of 
these same quantities occur naturally in conductivity equations. Although the units are 
not shown in equations (B8) and (B9), it is easily verified that H*i f? 0 H ii) and H*j 
(= VcPi]) are actually nondimensional. 

According to reference 7, equation (B7) can be expanded as 




V , T * 

x i x j H ij 


i=l H il 


TT* TJ* 

i=l j = l H ii H jj 


(B13) 


J* 


The first term in this expansion will be referred to as the first approximation 7]j and 
the first two terms will be referred to as the second approximation 77^. Thus 


v 2 

* _ y x f 
?i ■ Z H * 

i=l H ii 


(B14) 
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* 

% 


* x? 


V V TT * 

X.X.H4 


= y _ y y w 

A; -h* A/ A/ h* p]-* 

j=l 
lift 


Z—/ ^-7 TT* TT 3 * 

i=l "ii i=l j=l n ii n jj 


(B15) 


The translational thermal conductivity can be treated in a completely analogous 
fashion. According to Muckenfuss and Curtiss (ref. 13), the governing equation is 


X t = 4 


L 11 

L 12 

• ■ • 

L l, 

X 1 

L 12 

L 22 

• • • 

L 2j/ 

x 2 

L l, 

L 2i^ 

• • . 

h vv 

Xp 

X 1 

x 2 

. . . 

x v 

0 




(B16) 


where 


Lii = - 


X i X j 


* x i i6T y 

Xj 25p /-Li \2 

1=1 Mi + Mj) Dy 

jri ' ' 


[yMf -SB^ bS+^MjA*) 


(B17) 


16 t x iXj MiMj/gg A 

^ -5sF,.. . „\2 D„ (T 3B « 4A «) 


^ + Mjj U ij 


(i / j) (B18) 


In equation (B17), Xi, the thermal conductivity of pure component i, is defined by 


. J 75/ R’ 

1 = Wfe 


, 1/2 


X 


lol6 (ir) 2 


ll Wh 


'JJ'4 2 ’ 2T 


-(2,2) o 2 

where Oii is expressed in A . 


(B19) 
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After nondimensionalizing X(- with a reference conductivity defined by 

r ~ 


X Q - < 


75 B? 


a/2 


10 16 / 


mole 

lg-cm 4 


a/2 


Wt 




and using equation (B5) to eliminate Dy, equation (B16) can be expressed as 


where 


H*=r = 4 

A o 


T * 

L n 

T * 
^12 

. . . 


X 1 

L 12 

T * 
^22 

• • • 

L 2v 

x 2 

• 

• 

• 


* 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


T * 
L 2v 

• • • 

L* 

vv 

x ^ 

X 1 

x 2 

• • • 

x v 

0 




L *i = 


- ^ x i x j C 2 ,i j ^ 4 ,i j ^9 * 


5=1 

5/i 


(B20) 


(B21) 


(B22) 


C 4,ij ~ 


* -(1,1) 
L ij - x i x j c 2,ij c 5,ij n ij 

(i/5) 

(B23) 

f +4Ml MjA* 

Mj + Mj 


(B24) 


M i M j(f - 3B iJ ' 4A ij 

M. + M. 


Mj + 


(B25) 


With the assumed constancy of Ay and By, this formulation proves especially con- 
venient for numerical evaluation. 
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Equation (B21) can be expanded as 


v v 


x* - 
\ ~ 


V 9 V V _ * 

i=l Mi i=l j=l L ii L jj 

J* 


The first approximation is defined as 


* x? 




i=i L ii 


and the second approximation as 


v „2 


V V 


i=l Mi i=l j=l MiMj 

iA 


According to Butler and Brokaw (ref. 11), the reactive contribution to the 
tivity of a I'-component system whose equilibrium composition is described by 
pendent chemical reactions has the form 



0 

AH 1 

ah 2 


AH t 


AHj 

E 11 

E 12 


E lt 


ah 2 

E 12 

E 22 


E 2t 


• 

• 

* 


• 


• 

* 

• 


• 


• 

* 

* 


• 

1 

AH t 

E lt 

E 2t 


E tt 

RT 2 



E ik 




where AHj is the enthalpy of reaction of the ith independent reaction and 


(B26) 


(B27) 


(B 28 ) 


conduc- 
t inde- 


(B29) 
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V-l i-V / v \ 

V V * ( n kj x m “ n km x j) V V 2 A x r 

Ekk= 2 2 Xix m 2 2 "kj a Pi 3^ 

j=t+l m=j+l J j=t+l p=l J 

p^k 

y A ( x i+W 
A pkxk 1 x i xk 


2 v t 


i Xk • + i 

p=l ]=t+l 

P^k 


(k= 1,2,. ..,t) 
(B30) 


E V V A . (Hifm " n im x j)( n kj x m " n km x j) Y V n -- n , - A • ^ 

E ik - Z Z Am J . Z Z A PJ Xj 

i=t+l m=j+l ^ j=t+l p=l 


j=t+l m=j+l 

v 

+ 

j=t+l 


X ( n kj A ij + n ij A kj) " A ik 


a RT 
lj " pDij 


(i»k = l,2,...,t) 

(B31) 

(B32) 


Dividing equation (B29) by X Q and introducting equation (B5) for Djj yield 


0 

* 

AHj 

AH* 

• • * 

AHt 

* 

AHj 

E U 

E 12 

. . . 

E lt 

AH* 

E 12 

E ^2 

• • * 

E 2t 

AH* 

E *lt 

E 2t 

. . . 

E tt 


(B33) 


where 



(B34) 
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- 1 i 1 1 *** 

j=t+l m=j+l J j=t+l p=l 3 

pA 

+ y A . a + y a ..£jiw£ 

Lj Pk xjj Lj kj XjXt- 


2 v 


E 


kk 


P =1 

pA 


(k = l,2,...,t) 
(B35) 


v -1 v 


a. y y ~ -t mX K n kfm - Vi) f f 

/j /j m J XiX m Lj Lj PJ x-j 

3 m j=t+i p=i J 

i) - ^ 


j=t+l m=j+l 

v 

+ Z ( n kj A ij + “ij 4jj 

j=t+i 


A ii = — nff’^Cg u 

16 'MJ 


(i,k -= l,2,...,t) 
(B36) 
(B37) 


When equation (B33) is expanded, the first approximation is 



(B38) 
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TABLE 1 


BINARY INTERACTIONS FOR DISSOCIATING C0 2 -N 2 -Ar MIXTURES WITHOUT IONIZATION 


C0 2 -C0 2 C0 2 -0 2 

( 02 - 02 ) 


co 2 -n 2 

co 2 -no 

co 2 -co 

C0 2 -0 

g02-N| 

|c°2-c) 

C0 2 -Ar 

( 02 -N 2 ) 

o 2 -no 

o 2 -co 

M 

C°2-n) 

[02-c| 

0 2 -Ar 

(n 2 -n 2 ) 

(n 2 -no) 

n 2 -co 

(n 2 -o) 

(n 2 -n) 

n 2 -c 

N2-Ar 


(no-no) 

NO-CO 

(no-o) 

NO-N 

NO-C 

NO-Ar 



CO-CO 

co-o 

CO-N 

co-c 

CO-Ar 




1 

0 

(o-n) 

o-c 

O-Ar 





(n-n) 

N-C 

N-Ar 






c-c 

C-Ar 







Ar-Ar 


Parentheses around a pair designate that its interaction potential has been considered previously and that the 
corresponding collision integrals can be obtained from reference 19. 

Brackets around a pair mean that an accurate interaction potential is not critical because one or both of the 
mole fractions are always small. 
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TABLE 2 


INTERACTION POTENTIALS FOR DISSOCIATING C02-N 2 -0 2 -Ar MIXTURES 


Interaction 

Potential (r in A) 

Temperature range 

Reference 

Comments 

Interactions enclosed 
by parentheses in 
table I 

See reference 19 


19 

Averages of tabulated values were 
used for minor interactions 

CO-atom or molecule 

Use N 2 -atom or molecule potentials 



See section "Carbon monoxide 
interactions" 

Ar-0 2 

Exp-6: 0 = 15.5, e/k = 127.5° K, 
r xn - 3.796 A 

Low 

19, 25 

Obtained from viscosity data and 
combination rules 

1360 eV 

High 

29 

Scattering potential 


j.8.34 


At-N 2 

Exp-6: a = 15.5, e/k = 111.6° K, 
r m = 3.938 A 

Low 

25 

Obtained from viscosity data and 
combination rules 

76° eV 

High 

28 

Scattering potential 


r 7 * 78 



Ar-NO 

2382e~ 3 * 205r eV 

All 

18, 19, 20 

See section "Argon interactions" 

Ar-0 

1025e~ 3 ' 330r eV 

All 

18, 20 

Obtained from equations (10) 
and (11) by combination rule 


142 - 4 eV 

All 

28 

Obtained indirectly from scat- 


r 7.10 



tering data 

Ar-Ar 

Exp-6: a = 14, e/k = 123.2° K, 
r m = 3.866 A 

Low 

25 


32 300e” 4 - 46r eV 

High 

20 

Obtained from scattering data 

C02-C02 

6-12: 0 = 3.897 A, e/k = 213° K 

Low 

32 

Obtained from viscosity data 

678.4e" 2 * 166r eV 

High 

Appendix A 

Obtained by averaging 

CO2-O2 

Exp-6: a == 17, e/k = 255° K, 
r m = 3.616 A 

Low 

31 

Obtained from diffusion data 

320.8e' 2 * 153r eV 

High 

Appendix A 

Obtained by averaging; equa- 
tion (11) used for <£(0-0) 

C0 2 -N 2 

Exp-6: a » 17, e/k = 184° K, 
r m = 3.782 A 

Low 

31 

Obtained from diffusion data 

674.5e~ 2 - 4l6r eV 

High 

Appendix A 

Obtained by averaging; equa- 
tions (l4) and (15) used for 
<£(N-N) and <£(N-0) 


6-12: 0= 3.748 A, e/k = 139.2° K 

Low 

7 

Obtained from viscosity data and 
combination rules 

C0 2 -N0 

473.4e- 2289r eV 

High 

Appendix A 
...I 

Obtained by averaging; equa- 
tions (11), (14), and (15) used 
for <£(0-0), <£(N-N), and 
<£(N-0) 

C0 2 -Ar 

6-12: 0=3.681 A, eA - 157.2° K 

Low 

7 

Obtained from viscosity data and 
combination rules 

5401.7e~ 3 * 130r eV 

High 

Appendix A 

Obtained by averaging 

C0 2 -0 

151.3e~ 2 * 198r eV 

All 

Appendix A 

Obtained by averaging; equa- 
tion (11) used for <£(0-0) 

N-NO 

253.1e“ 2 ’ 647r eV 

All 

18, 19 

Obtained by averaging; equa- 
tions (14) and (15) used for 
<£(N-N) and <£(N-0) 

C-NO 

262.2e~ 2 ’ 670r eV 

All 

18 

Obtained by averaging; equa- 
tion (14) used for <£(N-C) 
and <£(0-C) 

C-0 

Use N-N potential 

All 


See section "Carbon interactions" 

Remaining 

C-Interactions 

Use N-interactions 



See sections "Carbon interactions" 
and "Carbon dioxide interactions” 


Potentials for CO?-N, CO2-C, O2-N, and O2-C are not included because either one or both of the species involved in these interactions 
are always small. Transport-property calculations were made with various potentials for these interactions and absolutely no change in the 
results occurred. 
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TABLE 3 


COLLISION INTEGRALS 


T , 

OK 


1000 

1500 

2000 

2500 

3000 

3500 

4000 

4500 

5000 

5500 

6000 

6500 

7000 

7500 

8000 


Ar 

-Ar 

Ar-C>2 

Ar-N 2 

Ar 

-O 

Ar 

-N 

Ar- 

•NO 

C02-0 2 

o M , 

a (2 - 2 >. 

a"- 1 ’, 

a< 2 ’ 2 >. 

a 11 ' 11 , 

n <2 " 2 >. 


S B,2) 


a< 2 ' 2 >. 

a* 1 ' 1 !, 


a* 1 ’ 1 *, 

a <2 ’ 2 >. 

A 2 

A 2 

A 2 

A 2 

A 2 

A 2 

A 2 

A 2 

A 2 

A 2 

A 2 

A 2 

A 2 

A 2 

8.73 

9.91 

8.94 

9.95 

9.40 

10.49 

6.98 

8.23 

6.90 

8.39 

9.03 

10.54 

9.60 

10.56 

7.50 

8.55 

8.22 

9.33 

8.56 

9.66 

6.35 

7.54 

6.16 

7.48 

8.29 

9.72 

8.94 

9.80 

6.92 

7.95 

7.61 

8.84 

7.80 

9.00 

5.93 

7.06 

5.67 

6.89 

7.78 

9.17 

8.50 

9.38 

6.48 

7.60 

7.13 

8.41 

7.18 

8.48 

5.61 

6.70 

5.33 

6.48 

7.41 

8.74 

8.20 

9.07 

6.20 

7.29 

6.78 

8.04 

6.75 

8.08 

5.36 

6.42 

5.08 

6.18 

7.10 

8.41 

7.98 

8.86 

6.04 

7.03 

6.54 

7.75 

6.48 

7.78 

5.15 

6.18 

4.83 

5.88 

6.85 

8.12 

7.80 

8.68 

5.88 

6.82 

6.35 

7.52 

6.27 

7.53 

4.97 

5.98 

4.67 

5.68 

6.63 

7.88 

7.65 

8.52 

5.76 

6.65 

6.17 

7.31 

6.08 

7.30 

4.82 

5.80 

4.51 

5.49 

6.45 

7.68 

7.51 

8.38 

5.65 

6.53 

6.02 

7.13 

5.92 

7.10 

4.68 

5.65 

4.38 

5.32 

6.28 

7.49 

7.38 

8.24 

5.55 

6.42 

5.88 

6.97 

5.78 

6.94 

4.56 

5.52 

4.27 

5.19 

6.14 

7.33 

7.26 

8.13 

5.46 

6.32 

5.76 

6.82 

5.65 

6.78 

4.45 

5.39 

4.16 

5.06 

6.01 

7.18 

7.14 

8.03 

5.38 

6.23 

5.64 

6.69 

5.54 

6.65 

4.35 

5.27 

4.06 

4.94 

5.89 

7.04 

7.02 

7.96 

5.31 

6.15 

5.55 

6.58 

! 5.44 

6.52 

4.26 

5.17 

3.99 

4.85 

5.78 

6.92 

6.91 

7.87 

5.23 

6.07 

5.45 

6.46 

! 5.34 

6.40 

4.17 

5.07 

3.92 

4.77 

5.68 

6.81 

6.80 

7.80 

5.17 

6.01 

5.36 

6.36 

| 5.25 

6.29 

4.09 

4.98 

3.83 

4.66 

5.58 

6.70 

6.70 

7.73 



co 2 - n 2 

co 2 -o 

C02-Ar 

co 2 -no 

co 2 - co 2 

NO 

-N 

NO-C 

T, 

°K 

A 2 

a< 2 ' 2 >, 

A 2 

A 2 

a l2 - 2 >, 

A 2 

a 11 - 1 *, 

A 2 

n (2 - 2 \ 

A 2 

A 1 

a (2 ’ 2 >, 

A 2 

A* 

a <2 ’ 2) , 

A 2 

a 11 ’ 1 ’, 

A 2 

a< 2 ’ 2 >, 

A 2 

A* 

n (2 - 2) 

°2 * 
A* 

1000 

9.94 

10.84 

9.72 

11.71 

10.90 

12.03 

11.05 

12.21 

13.08 

14.40 

7.84 

9.44 

7.78 

9.36 

1500 

9.26 

10.23 

8.75 

10.71 

10.14 

11.25 

10.30 

11.44 

12.10 

13.38 

7.02 

8.50 

6.96 

8.43 

2000 

8.84 

9.80 

8.01 

9.85 

9.68 

10.79 

9.88 

10.99 

11.50 

12.78 

6.46 

7.87 

6.41 

7.81 

2500 

8.54 

9.50 

7.46 

9.21 

9.28 

10.43 

9.51 

10.62 

11.13 

12.40 

6.05 

7.39 

6.00 

7.34 

3000 

8.30 

9.28 

7.02 

8.70 

8.93 

10.15 

9.14 

10.31 

10.83 

12.05 

5.72 

7.01 

5.68 

6.97 

3500 

8.10 

9.08 

6.65 

8.28 

8.64 

9.89 

8.78 

10.07 

10.55 

11.78 

5.45 

6.70 | 

5.41 

6.66 

4000 

7.94 

8.90 

6.35 

7.93 

8.39 

9.68 

8.46 

9.86 

10.30 

11.52 

5.22 

6.44 

5.19 

6.40 

4500 

7.79 

8.76 

6.09 

7.62 

8.18 

9.51 ; 

8.18 

9.67 

10.04 

11.30 

5.02 

6.21 

4.99 

6.17 

5000 

7.65 

8.63 

5.86 

7.35 

7.99 

9.35 

7.93 

9.48 

9.80 

11.15 

4.85 

6.01 

4.82 

5.97 

5500 

7.51 

8.51 

5.67 

7.11 

7.82 

9.20 

7.69 

9.30 

9.56 

10.98 

4.69 

5.83 

4.67 

5.79 

6000 

7.37 

8.41 

5.46 

6.90 

7.67 

9.05 

7.49 

9.11 

9.34 

10.88 

4.55 

5.67 

4.53 

5.63 

6500 

7.24 : 

8.31 

5.31 

6.70 

7.53 

8.90 

7.30 

8.94 

9.12 

10.75 

4.43 

5.52 

4.40 

5.49 

7000 

7.11 

8.22 

5.16 

6.52 

7.40 

8.76 

7.13 

8.77 

8.92 

10.66 

4.31 

5.39 

4.29 

5.36 

7500 

6.98 

8.16 

5.02 

6.36 

7.29 

8.64 

6.97 

8.61 

8.73 

10.53 

4.21 

5.26 

4.19 

5.23 

8000 

6.85 ! 

8.08 

4.89 

6.21 

7.17 

8.51 

6.83 

8.44 

8.56 

10.43 

4.11 

5.15 

4.09 

5.12 
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(a) p = 0.001 atm. 

Figure 1.- Equilibrium mole fractions of high-temperature air for several pressures. 
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Species mole fraction 
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Species mole fraction 



Temperature, °K 

(c) p = 0.1 atm. 
Figure 1.- Continued. 
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I 



(d) p = 1 atm. 
Figure 1.- Continued. 
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Temperature, °K 

(e) p = 10 atm. 
Figure 1.- Concluded. 
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lb) p = 0.01 atm. 
Figure 3.- Continued. 
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(c) p = 0.1 atm. 
Figure 3.- Continued. 
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(d) p = 1 atm. 
Figure 3.- Continued. 
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Temperature, °K 


(e) p = 10 atm. 
Figure 3.- Concluded. 
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Figure 4.- Variation of nondimensional mixture specific heat of air with temperature. 
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Figure 6.- Variation of nondimensional mixture specific heat of Mars model atmosphere with temperature. 
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5.0 



Temperature, °K 
(a) O2 and N2. 

Figure 8 .- Variation of nondimensional species enthalpy with temperature. 
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T empe rat ur e , 0 K 


(c) C, N, 0, NO, and Ar. 
Figure 8.- Concluded. 
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Figure 10.- Comparison of Q* 1 ’ 1 * 


computed from high-temperature Ar-Ng and Ar-CO potentials. 
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: igure 17.- Illustration of graphical method of finding 


n <u> 


for Ar-02 interaction. 
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Figure 18.- 


I llu stration of graphical method of finding 


—(2 2 ) 

Q ' for Ar-02 interaction. 


> 0 K 











Figure 24.- Variation of nondimensional total conductivity of air with temperature. 


89 




1 2 3 4 5 6 7 8x 10- 

Temperature, °K 

Figure 26.- Comparison of full-equation calculation of nondimensional viscosity of air with approximate calculations 

and previous results for p = 1 atm. 




Figure 27.- Comparison of full-equation calculation of nondimensional frozen conductivity of air with approximate calculations 

and previous results for p = 1 atm. 
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Figure 28.- Comparison of full-equation calculation of nondimensional reactive conductivity of air with approximate calculations 

and previous results for p = 1 atm. 




Figure 29.- Comparison of full-equation calculation of nondimensional total conductivity of air with approximate calculations 

and previous results for p = 1 atm. 






Temperature, °K 

Figure 31.- Comparison of full-equation and approximate calculations of Prandtl number of air with modified data of Yos for p = 1 atm. 
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Figure 34.- Variation of nondimensional reactive conductivity of CO2 with temperature. 
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Figure 37.- Comparison of full-equation calculation of nondimensional viscosity of CO 2 with approximate calculations 

and previous results for p = 0.1 atm. 
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Temperature, °K 

Figure 38.- Comparison of full-equation calculation of nondimensional frozen conductivity of CO 2 with approximate calculations 

and previous results for p = 0.1 atm. 
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Figure 39.- Comparison of full-equation calculation of nondimensional reactive conductivity of CO 2 with approximate calculations 

and previous results for p = 0.1 atm. 




Figure 40.- Comparison of full-equation calculation of nondimensional total conductivity of CO 2 with approximate calculations 

and previous results for p = 0.1 atm. 
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Figure 42.- Variation of nondimensional viscosity of the Mars model atmosphere with temperature. 
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Figure 44.- Variation of nondimensional reactive conductivity of the Mars model atmosphere with temperature. 
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Figure 46.- Variation of Prandtl number of the 
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Figure 47.- Comparison of nondimensional viscosities of air, CO 2 , and the Mars model atmosphere for p = 1.0 atm. 
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Figure 48.- Comparison of nondimensional total conductivities of air, CO 2 , and the Mars model atmosphere for p = 1.0 atm. 
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Figure 49.- Comparison of Prandtl numbers of air, CO 2 , and the Mars model atmosphere for p = 1.0 atm. 
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